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Abstract 

The progress of image processing during recent years allows the measurement of 
pedestrian characteristics on a "microscopic" scale with low costs. However, density 
and flow are concepts of fluid mechanics defined for the limit of infinitely many 
particles. Standard methods of measuring these quantities locally (e.g. counting 
heads within a rectangle) suffer from large data scatter. The remedy of averaging 
over large spaces or long times reduces the possible resolution and inhibits the gain 
obtained by the new technologies. 

In this contribution we introduce a concept for measuring microscopic character- 
istics on the basis of pedestrian trajectories. Assigning a personal space to every 
pedestrian via a Voronoi diagram reduces the density scatter. Similarly, calculating 
direction and speed from position differences between times with identical phases of 
movement gives low-scatter sequences for speed and direction. Closing we discuss 
the methods to obtain reliable values for derived quantities and new possibilities 
of in depth analysis of experiments. The resolution obtained indicates the limits of 
stationary state theory. 

Key words: Video tracking, Voronoi diagram, pedestrian modeling, velocity 
measurement, pedestrian density 



1 Introduction 



For the design of pedestrian facilities concerning safety and level of service 
piT21[5Ello] . pedestrian streams are usually characterized by basic quantities 
like density or flow borrowed from fluid dynamics. Up to now the experi- 
mental data base is insufficient and contradictory [3|6|7] and thus asks for 
additional experimental studies as well as improved measurement methods. 
Most experimental studies of pedestrian dynamics use the classical definition 
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of the density in an area by D = r^r, where N gives the number of pedestrians 
in the area A of size \A\ [m 2 ]; see e.g. [8ll9lfT0lfTTir^^ 

An obvious problem is that this gives an averaged density for a specific area, 
not a density distribution p(x) [1/m 2 ], and passing to the limit of areas of 
size zero is obviously not a well defined procedure. Also the choice of the 
geometry of A is important. For large convex areas it can be expected that 
finite size and boundary effects as well as influences of the shape of the area 
can be neglected, though it is almost always possible to cut out fairly large 
areas (of complicated shape) containing no person. Further, the design of 
pedestrian facilities is usually restricted to an order of magnitude of 1 to 10m 
and rectangular geometries. In addition the number of pedestrians N is small 
and thus the scatter of local measurements has the same order of magnitude 
as the quantity itself, see e.g. the time development of the density in front of 
a bottleneck in [15] . Often, due to cost restrictions the density was measured 
at a certain point in time, while the measurement of speed was averaged over 
a certain time interval [9XTT] . But the process of measurement and averaging 
influences the resulting data even for systems with few degrees of freedom as 
the movement along a line, see [T9.22J. 

The large progress in video techniques during recent years has made feasible 
the gathering of much more detailed data on pedestrian behavior, both in 
experiments and in real life situations, than was possible only a decade ago 
[T5|20|21j . The higher detail asks for a reevaluation of the methods of defining 
and measuring basic quantities like density, flow and speed, as the methods 
to get time and space averages encompassing a hundred persons over minutes 
may not be suitable for a resolution of a second and a square meter. Basic 
quantities of pedestrian dynamics are the density D [1/m 2 ] in an area A, the 
the velocity v [m/s] and speed s = \\v\\ of a person or a group of persons, 
and the flow through a door or across a specific line F [1/s] . Measurements 
give averages of these quantities, and trying to carry the measurement to the 
infinitesimal limit is obviously not reasonable. Below we indicate where the 
limits for the possible resolution are, and give methods that allow to go fairly 
close to these limits. The methods presented here are based on video tracking 
of the head from above, but tracking of e.g. a shoulder or the chest might do 
even better, though they are more difficult to obtain. 

As density is used as one (important) parameter in modeling peoples behavior, 
a perceived density would be best, but the details of perception are completely 
unknown. However, the perception definitely does not show the rapid fluctu- 
ations in time that the classical measurements do. 

There is not a 'right' definition for density, flow or velocity of pedestrians 
(though some definitions may be wrong), the definitions treated here are more 
or less useful for a certain purpose and allow more or less resolution. 
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All the examples in the paper are from experiments conducted in Jiilich and 
Diisseldorf and described in [21] and [22J. 



2 Direct Measurements 

2.1 Measuring density 

In pedestrian dynamics, density is persons per area. The density at a point is a 
mathematical abstraction, as persons are discrete and extended in space, and 
the density is well defined only on a scale large enough to be able to ignore 
discreteness and small enough to be homogeneous. These two conditions are 
no problem in fluid mechanics with > 10 18 particles per mm 3 , but in treating 
pedestrians they are conflicting. The definition is directly reflected by the 
standard procedure of measuring it. Here a certain area, usually a rectangle, 
is laid out, and the number of persons inside this area is counted. With N 
pedestrians in the measurement area A the standard definition is 



This method has two drawbacks. The less important one is that occasionally 
"in" and "out" has to be assigned arbitrarily, usually by head position. More 
important, the density depends discontinuously on time and on the exact 
placement of the measuring area, with large jumps for small area. These jumps 
can be reduced by taking averages over time and/or position of the area, at 
the cost of resolution. 

A way to define point values of the density is having every person i pro- 
duce a density distribution Pi(x) - a non-negative function with unit integral 
/ Pi(x) dx = 1. A step function (e.g. p%{x) = l/(2r 2 7r) for \\x — xA\ < r) , a 
linear function of the distance Pi(x) = max(0, h(r — \\x — Xi\\), or a Gaussian 
may do for different purposes, e.g. in molecular dynamics. The central prob- 
lem is the extention r (or a for Gaussian) of this function. The requirements 
are that in a situation considered homogeneous, the resulting density from a 
group of people will not show too much variations, while inhomogeneities will 
not be masked by too large a width of the individual density function. These 
requirements can't be met with a fixed extention r (or a) of Pi. Much more 
involved but no problem for computers, is the following adaptive procedure 
for step functions: 

• With a given set of trajectories of M persons {x\(t), x^it), XM(t)} in two 
dimensions x = (x,y) assign an exact position Xi(t Q ) at time t for every 
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person i. 

• Compute the Voronoi-diagram [23] for these positions, giving cells A± for 
each person i. 

• Compute the size of the cell \Ai\ = J A dx and define the density distribution 
for all persons: 



Pi(x) 



W\ : x e Ai 
: otherwise 



and p(x) = ^Pi 



x 



(2) 



The Voronoi diagram FigJT] assigns to every point the area closer to this point 
than to any other one. It has been used for defining density in the fields of 
granular media [24], molecular physics [23] and neuroanatomy [25], and finds 
further widespread applications in graphics and mesh generation for PDEs. 
Here the program VRONI is used [27]. 



x 

Fig. 1. Positions of persons jamming in front of an exit with exemplary Voronoi 
cells 



Problems arise for the persons at the rim of a group, for which the Voronoi 
cell may extend to infinity. Restricting the individual cells in size helps, but 
the details of this are somewhat arbitrary. For our procedure to measure the 
density in the given area, we used a restriction to 2m 2 , which is active only for 
a few cells. In the presence of walls, the Voronoi cell is restricted to the area 
inside. It has been proposed by Pauls [28] for stairs to use virtual walls 15 cm 
inside the physical walls, but this distance depends on details of the walls like 
roughness etc. and may be different for level floors. 

With the procedure defined above different density definitions are possible. To 
define a density in an area A we choose 

„ [ A p(x) dx , , 

Dv = JAP ^ (3) 
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Another Voronoi diagram based definition which has some merits with respect 
to the observation area required is 



N 



Dy = _ N n , (4) 



In this definition only the sizes of the Voronoi cells for the persons inside A 
are used. Both Dy and Dy are not very sensitive to the positioning of the 
area of measurement. Dy is preferable for small measurement areas inside a 
crowd, where Dyi may not be defined if there is no person inside the area. 
On the other hand, Dy requires useful Voronoi cells for all people whose cell 
overlaps with the area A, not just for those whose head is inside A as Dy 
does. Therefore it is more likely to be affected by proximity to the rim of a 
group. 




»[m] 



(a) Classical density D s (b) Voronoi density Dy (c) Voronoi density Dy 

Fig. 2. Different definitions of densities in an observation area A. Definition of D s is 
according Eq. [fl while Dy and Dy is according Eq. [3] and Eq. H] respectively. The 
dark blue area contributes to the densities Dy and Dy. 

Fig. Ogives a time sequence of D s and Dy. The latter is much more stable in 
time, with less than half the standard deviation. While D s changes repeatedly 
by about 30% and back within a second, there is only one fast change in Dy, 
and this is less than 10%. Thus while D s at a single time gives no useful 
information, Dy does. 

Fig. [3] shows that the overall picture is similar for all definitions of densities. 
All densities are measured in the square A of \A\ = 1 m 2 , see Fig. [2j D s 
obviously varies much more than the others, and more than the perception of 
the situation of a human observer, because entering or exiting of the rectangle 
by a person changes D s considerably, but not the perception. (D s ) space was 
averaged over rectangles shifted by up to 50cm in each direction, {D s ) time is a 
moving average over 21 frames (0.84s). (D s ) time shows more variation in time 
than (D s ) space , which again varies more than Dy. Dy has the advantage of 
having full resolution in space and time combined with low fluctuations. 
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Fig. 3. Time sequence of standard density D s and Voronoi density Dy. (a): Total 
run; (b): 'Quasi' stationary state t = [20,32]s; (c): D s , {D s ) space and (D s ) t i me ; (d): 
D s , By and (D s ) time . 



For the different densities from Fig. 3 we have the data in Table [TJ The differ- 
ence between the density averages for D s and Dy are within the limits of the 
fluctuations, see Table [fl but it is obvious that the density distribution is not 
homogeneous over the entire camera area. The Voronoi cells carry information 
from outside the rectangle, where the density may be different. This may be 
the reason for part of the differences. 

Using Voronoi cells, a density distribution p(x) is attributed to every point in 
space. However, this distribution oscillates with stepping, so for best results 
only time averages have to be taken over the time of at least a step, or some 
smoothing of the oscillations is needed (see below). 
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Table 1 

Mean value O, standard deviation c(0) and total variation TV{0) = ^ ■ |Oj+i- 
for different density definitions for results from Fig. [31 t = [20, 32] s. 

2.2 Measuring velocity 



Given trajectory Xi(t) of a person i, one standard definitions of the velocity 
with At fixed, but arbitrary is 



VAt,i(t) 



Xi{t + At/2) -Xi(t- At/2) 
At 



(5) 



Alternatively, with given entrance and exit time ti n ,t out the velocity is 



tout tin 



(6) 



with s = \\v\\. The average standard speed in an observation area A is then 

s s (t) = l/N s At,M (?) 

Xi(t)eA 



where the sum is taken over all persons that are in A for the entire time in- 
terval [t — At/2, t + At/ 2]. These definitions seem simple, but there are two 
sources of uncertainty. The velocity of an extended object is generally (and 
reasonably) defined as the velocity of its center of mass, and for pedestrians 
that is hard to detect. Moreover, pedestrians can and do change shape while 
walking. So the simple approach works only for distances long enough to make 
errors from shape changes and in placing the supposed center of mass unim- 
portant. A second problem comes from the fact that velocity is a vector, and 
the movement of people is not straight. Thus the average of the local speeds 
will be bigger than the value of distance per time for longer distances. No- 
tably head tracking gives tracks that can be decomposed into a fairly uniform 
principal movement and a local swaying superimposed. The swaying shows 
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the steps, it varies between individuals and is larger at low speeds. Of course, 
the head sways more than the center of mass, and it can do some independent 
movements, but these can dominate only at very low speeds. For most mod- 
els of pedestrian movement, only the principal movement is of interest. The 
swaying movement of shoulders may be important in estimating the necessary 
width for staircases and corridors to allow overtaking. 

The separation of principal movement from swaying movement could be done 
by Fourier analysis, but that requires a trajectory many steps long. A way 
to do it locally is to detect positions of identical phase of the movement and 
interpolate them. As long as there is appreciable forward motion (> 0.3 m/s), 
the mode of movement is the swinging of the legs in the direction of move- 
ment with approximately their free pendulum frequency (1.5Hz-2Hz). In this 
mode there is a regular sequence of points of maximum (positive) curvature, 
minimum (negative) curvature, and zero curvature, which correspond to the 
times of setting down the right foot, the left foot, and having one foot on the 
ground while the other just passes the standing leg in forward motion, which 
are the points we use. These points are easy to detect. Below that speed, the 
mode of stepping changes to the whole body swinging right and left with a 
frequency smaller than 1Hz and only a small forward component, and there 
may be multiple points with zero curvature d 2 Xi = within one step. How- 
ever, typically the positions of setting down a foot give a dominant extremum 
of curvature, while in the part in between the curvature will be close to zero 
with more than one zero per step. In this case, we take the middle between 
the maximum and the minimum curvature point as interpolation point. This 
has been possible down to speeds of 5 cm/s. Below that, steps can only be 
guessed, they can not be detected reliably. 

The speed of the principle movement s Pj j = \\v p ,i\\ is calculated by 



Interpolating these points now eliminates most of the swaying movement and 
gives a good approximation to the movement of the center of mass. The re- 
quirement of identical phase asks for taking only every other zero curvature 
position, but for persons with symmetric gait, taking every zero curvature 
position gives a better resolution with only marginally more swaying. For 
analysis, we will take this as curve of principal movement. 

Similarly, the velocity vector will be obtained by computing the difference 
quotient of position and time between zero curvature positions, and attached 
to the intermediate time. 



v, 




with curvature(xj(^)) = 0, j — 1,2. (8) 



*2-*l 
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x [m] 

Fig. 4. Trajectories of three persons walking through bottleneck. Circles denote 
interpolation points for a smoothed trajectory. 
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Fig. 5. Speed in x and y versus time for part of the middle trajectory of Fig. 
dotted(red): At = 0.04s; dashed(black): At = 0.44s according Eq. [5j solid(blue): as 
described in Eq. El 

2. 3 Measuring flow 



Standard measurement of flow F s through a door or across a line is done similar 
to density measurements, by counting heads passing within a time interval. 
This suffers from the same problems as the standard density measurement, 
namely large scatter and low time resolution. Using the Voronoi cells to obtain 
fractional counts (half a person has passed if half of the Voronoi cell has 
passed) gives a much smoother Voronoi flow Fy. This still does not allow a 
useful passing to the limit of At = 0, but the moving average over about half 
the average time difference between persons passing gives a sufficiently smooth 
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result. 
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Fig. 6. Flux through bottleneck with standard method F s and Voronoi cell based 
method F v , left: frame by frame right: 1 s moving average. 
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Fig. 7. Flux at bottleneck. Voronoi method F v and flow computed by F = D s s s b 
where b is the width of the bottleneck. 

Fig.Hcompares the product of averaged standard speed and density, s s (t) D s (t), 
with the direct measurement of the flow as Voronoi density passing the middle 
line. The density was determined from a rectangle symmetric to the line for 
the flow of width Ax and the averaged speed of the persons inside this rect- 
angle from a symmetric time difference of At. The two product evaluations 
agree reasonably well with the Voronoi flow, but in spite of the fact that the 
measurements average over time and space they show much faster variations 
in time, and the total flow they calculate is somewhat less than the number of 
persons passing. For Ax = lm, At = 0.32s one person of 180 is missing in the 
integration, for Ax = 2m, At = Is it is even 7 persons missing. Depending 
on the purpose of the measurements, this may be a serious problem. 
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3 Prospects of analysis 

3.1 Fundamental diagram at bottlenecks 




The high resolution measurement of speed and density allows to follow indi- 
viduals on their way through some obstacles and look what combinations of 
speed and density they have. In Fig. [8] the persons start in front of the bottle- 
neck with low speed, no. 72 and 80 with large space while no. 58 is already in 
the jamming area. They pass the congestion area and pick up speed inside the 
bottleneck. It also allows correlating momentary speed and personal space, as 
well as comparing personal space and speed in general for whole groups of 
persons. Fig. [8] and [9] show a substantial difference in this relation before and 
inside the bottleneck, indicating that the individual speed depends more on 
the expectation of the future (walking into or out of high density regions) than 
on the present situation. The high resolution makes it also possible to analyze 
how the accelerations observed are related to changes in the space available, 
and to the space of people in front. On the resolution presented here it be- 
comes clear that the correlation of speed and available space in an instationary 
situation differs considerably from that in the stationary situation described 
by the fundamental diagram. 



11 




The most important relation for any model (and for much of the analysis 
of pedestrian movement) is the so-called fundamental diagram, which can be 
given either as relation speed versus density or flow versus density. Using the 
improved methods of measurement, the quality of the resulting diagram is 
greatly enhanced. 
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Fig. 10. Speed versus density for analysis with different methods. The trajectories 
are from an experiment of single file movement, see [22]. Left: sampled when a 
person crosses a certain line. Right: sampled every second, high densities only. Di- 
amonds(red): SAx,i with Ax = 2m versus density D s . Crosses (blue): same speed, 
density Dy. Points (black): speed from step detection s Pt i, density Dy. 

As an example Fig.f fTOl) shows the analysis of a series of measurements of 
the fundamental diagram for single file movement performed with different 
numbers of persons (14 to 70) in the walking area, which covers all densities 
of interest from low (almost free walking) to jamming density (near standstill). 
With the high density, stop-and-go waves developed, see Fig. 5 in [22j . For the 
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one dimensional case with position Xi of pedestrian i the calculation of the 
Voronoi density distribution reduces to 



Pi(x) 



&i — 1 



x e [ 



otherwise 



3/j — \~r^i Xj i J r \-\-X<i 



(9) 



The definition of density Dy is according to Eq. [31 Fig. (flU!) (left) shows 
that use of Dy greatly enhances the quality of the diagram, due to the fine 
density resolution. Use of the standard density contracts the band of s versus 
D onto a few vertical lines which are much longer than the width of this band 
and thus reduces the precision. In Fig. ffTO]) . (right) one can see that with the 
standard speed SAi,i the lowest speed is about 0.02 m/s, while s P)i gives values 
of zero. Actually, some people were standing in one position for about 30s - the 
duration of a stop phase. For higher speeds, there is little difference between 
$Ax,i and s Pj i. The combination of Dy and s Pt i gives the best diagram for the 
full scale of densities. 



4 Concluding remarks 



The combination of modern video equipment with new methods for extracting 
relevant data allows an unprecedented depth of analysis of pedestrian behav- 
ior. The method for determining density is based on the concept of a Voronoi 
cell as personal space of a pedestrian and allows a resolution down to indi- 
vidual level. The concept of determining velocities from difference quotients 
of positions with identical phase of the stepwise movement gets the resolution 
down to a single step. This level of resolution allows mathematical combina- 
tion of data that are not valid for large scale averages. Moreover they are able 
to resolve stop and go waves and allow a analysis of instationary processes on 
a microscopic level. 
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